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Abstract--The electromagnetic fields of point sources with time varying charges moving in the vacuum 
are derived using the Liénard-Wiechert potentials. The properties of the propagation velocities and the 
Doppler effect are discussed based on their far fields. The results show that the velocity of the electromagnetic 
waves and the velocity of the sources cannot be added like vectors; the velocity of electromagnetic waves of 
moving sources are anisotropic in the vacuum; the transverse Doppler shift is intrinsically included in the 
fields of the moving sources and is not a pure relativity effect caused by time dilation. Since the fields are 
rigorous solutions of the Maxwell’s equations, the findings can help us to abort the long-standing 
misinterpretations concerning about the classic mechanics and the classic electromagnetic theory. Although 
it may violate the theory of the special relativity, we show mathematically that, when the sources move faster 
than the light in the vacuum, the electromagnetic barriers and the electromagnetic shock waves can be clearly 
predicted using the exact solutions. Since they cannot be detected by observers in the region outside their 
shock wave zones, an intuitive and reasonable hypothesis can be made that the superluminal sources may be 
considered as a kind of electromagnetic blackholes. 


I. Introduction 


Maxwell’s theory is the foundation for handling all electromagnetic problems. In the coordinate 
system O (r,t), the Maxwell’s equations in the vacuum can be expressed by [1-3] 


V xE(r,7) =~ B(r,1) 
VxH(r,1)= PED 40 o 
V-B(r,t)=0 


V-D(r,t)= p(r,t) 


where E(r,ż) is the electric field intensity, H(r,t) is the magnetic field intensity, D (r,t) = £,E (r,t) 
and B(r,t) = 4H (r,t) are respectively the corresponding flux densities. £, and yz, are respectively the 
permittivity and the permeability of the vacuum. The fields are generated by the charge density p(r,t) 
and the related current density J (r,t). In the vacuum, the current density is caused by the motion of the 
charge density, i.e., 


J(r,¢) = p(r,t) V(r, 6) (2) 


v(r,t) is the velocity of (r,t). In this paper, we consider non-relativistic electromagnetic problems, 


and ignore the constraints from the theory of special relativity (SR). In the analysis, we are supposed to 
stay in the same coordinate system O (r,t) and do not care about the observers in other inertial frames. 


In the Maxwell’s equations, ¢, and 4, are constants in the vacuum, hence, c, = 1/ af L£, is also a 


constant. However, it is obvious that the Maxwell’s theory does not state that the velocity of the 
electromagnetic waves in the vacuum is always c,. Therefore, we temporarily forget the relativity rule 


that all objects cannot move faster than the light in the vacuum, and consider that all causal solutions to 
the Maxwell’s equations are physically reasonable. 


The electromagnetic fields generated by point charges in the vacuum are very useful because they 
usually have explicit expressions that can be used to illustrate the main characteristics of the 
electromagnetic fields. Liénard-Wiechert potentials [1, 4, 5] are widely used in analyzing the 
electromagnetic fields of a moving charge. They are derived by Liénard in 1898 and Wiechert in 1900. 
The same techniques can be applied for deriving the fields generated by moving sources with time- 
varying charges, like the Hertzian dipole [3, 6, 7, 8]. However, most of the published results for moving 


Hertzian dipoles are relativistic ones that were used to check the relativistic behaviors. In this paper, the 
rigorous solutions of a moving point source with time-varying charge and that of a moving Hertzian 
dipole are derived based on the Liénard-Wiechert potentials. The behaviors of the far fields are analyzed 
and illustrated with figures. In particular, the exact relationship between the wave velocity and the 
velocity of the sources are provided, based on which the anisotropic property of the wave velocity and 
the Doppler effect are demonstrated. These findings are very important because they have broken the 
three long-standing misinterpretations concerning about the classic physics. Statements can be found in 
most of the related textbooks and Journal papers that, according to the classic physics, the light velocity 
is isotropic in the vacuum; the wave velocity and the velocity of the source can be added like vectors; 
the transverse Doppler shift is a pure relativistic effect due to the time dilation. We show that all these 
statements are not true. Furthermore, without the velocity limit imposed by SR, when the sources move 
faster than the light in the vacuum, the electromagnetic barriers and electromagnetic shock waves can be 
clearly predicted from the exact solutions. 


II. Fields of Moving Point Sources 


Assume that a particle with harmonic charge COS@j, moves along a trajectory x(t) , as shown 
in Fig. 1. @ is the oscillating angular frequency of the charge. At the time ż,, the velocity is 
v(¢,) =x(t,), and the acceleration is a(t,) = v(t, ) . For the sake of brevity, we use the dot on top of a 
vector to denote the derivation with respect to time t,. Denote R =r —x(ż, ) as the radius vector from 
the source position to the observation point, and R = |r -x(t ) . The fields generated by the charge at 
time ¢, propagate to the observation point r with a time delay of t—¢,. In this paper, we generally use 
(r,t) for the fields and potentials, and (r,,t,) for the sources if not specified otherwise. For the sake of 


convenience, we assume that the point source at time ¢ is near to but not fixed at the origin, so that we 
can reveal the temporal behavior of the far fields. 


trajectory 


R 


E(r,t) 
Fig.1. Trajectory of a moving point source with time-varying charge. 


The charge density and the corresponding current density can be expressed by 
P(r,.t,) = 6 (1, —x)cos(@t, ) (3) 
I(r,,t,) = y(r, —x)cos( apt, ) (4) 


where (t,) in v(t,) and x(t,) is suspended for the sake of simplicity. Following exactly the same way 


in deriving the Liénard-Wiechert potentials for a moving charge [4, 5], we can derive the Liénard- 
Wiechert potentials for the moving point source with harmonic time varying charges, 


2 oo 
Drar (r,t) = ATE, R(1-f-B) cos(@t,) (5) 
A (r,t) = Po B cos (@t, ) (6) 


47e cyR(1-i-B) 


From the potentials, we can derive the fields of the moving source as [4, 5] 


= OA par (r, t) 


E (t,t) = a -VØ (r,t) 
ta, pe A cos( œt ) (7) 
_ Py R’ (1-n-B) c,R(1-n-B 
freo | AB ~ @, Sin (ot) 
cR(1-n-B 
Bua (r,t) = Lâ X Ear (r,t) (8) 
c 


0 


In the expressions, we denote B = v/c, , È = v/c, À =R/R, and c, is the light velocity in the vacuum. 
According to the solutions of the Maxwell’s equations, the field at (r,t) generated by the pulse source 


at x(¢,) travels with velocity c,. Hence, we have 


R=c,(t-t,)=|r—-x(¢,)| (9) 


When a, =0, it is straightforward to check that (7) and (8) are simplified to the fields of the moving 
point source with a constant charge. 


When the source moves uniformly, x(¢,)=v/, and B(¢,)=0. Denote B = BB, r=rt. B and ĉ 


are the corresponding unit vectors, and B-f =cos@. Solving ¢t, from (9) yields 


Cot —¢, Br cos 8 -c di-e sin? 0) r? —2c,tBcosOr + B’cct? 
O co (1- 8°) 


Note that ¢, must be real, and ¢, <¢ must hold according to the causality principle. For a uniformly 


(10) 


moving source, we can check that there is at most one solution for ż, that satisfies these conditions if 
2 <1. In other words, the fields at the position r and the time instant ¢ comes only from the source at 
the time instant ¢, when it moves to the position wt,. We may consider the field generated at a point 
(r,,¢,) as a pulse of electromagnetic fields. The fields generated at each time instant ¢, will propagate 
with constant velocity c, . However, a continuously moving point source generates continuous 


electromagnetic fields when it moves on. All the field pulses by the moving point source are superposed 
to form the continuous field distributions in the vacuum. However, we can see that the behaviors of the 
composed fields may become quite different from that of a single electromagnetic pulse: both the wave 
velocity and the frequency of the composed fields will change. 


At places far away from the sources, r >> vt , we have 


c, cosé 


J1- 8’ sin’ 0 


d-e sin? 0) r? —2c,tBcosOr + Bot” ~ ae sin’ Or — 


Consequently, 


(11) 


j l 1 B cos 0 } BcosO+J1— B sin’ 0 
1 > r 


+ 2 
1- f° 1-8 Ji- 6? sin? @ c(1- 8°) 


Substituting (11) into (7) and discarding the term with respect to (1/ R°), we obtain the approximate 


expression for the electric far field 


a Po n-B H = 
E,,, (r,t) * eae Rap sin| a(t r/c) | (12) 


where the velocity c and the angular frequency @ of the composed far field are, respectively, 


es a (13) 
1- 8’ sin’ 8 
wz O j 2s (14) 
1-8 1- 2’ sin’ 0 


Obviously, (12) represents a spherical wave. The wave velocity c and the angular frequency @ of the 
composed far fields are both functions of three variables: the velocity 2 of the source, the oscillating 


frequency @, of the source, and the angle 0 between the propagation direction of the fields and the 


moving direction of the source. In particular, equations (12)-(14) are valid for fields at regions far away 
from the point source. They do not explicitly dependent on the coordinate of the point source, therefore, 
it is convenient for us now to put the point source at the origin and analyze its far fields with (12)-(14). 


A. Wave velocity 


The wave velocity of the electromagnetic far fields is dependent on the propagation direction. It is 


obvious that c> c, for all 0. The velocity reaches its maximum of c, / Vl-f# at @=+7/2. It is 
exactly c, at @=0 or z, in which the observer is on the path of the source. The wave velocity of the 
far field is isotropic only in the case that 2 = 0, i.e., the source is motionless or moves within a bounded 
region. To state that the wave velocity in the vacuum is constant according to the classic physics rules is 
a misinterpretation. 

A typical plot of the velocity is shown in Fig.2. When £ = 0.1, the largest wave velocity is only 
1.005c, by (13). It would be 1.1c, if the two velocities are added directly. Equation (13) clearly shows 


that the propagation velocity of the electromagnetic waves and the moving velocity of the sources do not 
satisfy the vector addition rule. To state that the two velocities are added like vectors according to the 
classic physics rules is another misinterpretation. 


41-8 


Fig. 2. The velocity of a uniformly moving source. 


B. Doppler effect 


The angular frequency @ of the far fields is also dependent on the propagation direction. This is the 
Doppler effect. We denote the normalized angular frequency as 


1 B cos 0 


s(0)=—=—+ (15) 
o, l-8 1-6 Ji- p sin @ 
The relative Doppler shift is s(@)—1. The Doppler shift is zero at the angle @, that satisfies 
cos@, =- p (16) 


V1+ 8? 
For slowly moving sources, is small. s (0) can be expanded to the power series with respect to £ , 
8(0) ¥1+ BcosO+ B? +- (17) 


We can check that s(0)=1+ 6/(1- 8) 1+ and s(z)=1+ B/(1+ B)*1-f. They agree with the 


classic mechanics results and the results based on SR to the first order of £. 


Particularly, (14) shows that transverse Doppler effect clearly exists at 9=+7/2, 


1 2 

The Doppler effect is explicitly included in the wave solutions of the Maxwell’s equations at places 
far away from the sources. It is caused by the superposition of the fields radiated by the moving source 
at different times. It is directly solved in the laboratory frame without necessity to apply any invariant 
principles related to coordinate transformations. We have to note that the coefficient of the second-order 
Doppler shift is different from that predicted using SR, which is (—0.5) [4]. 


Conventionally, the transverse doppler effect is considered as a pure relativistic effect that is directly 
related to time dilation [4, 9]. It has been used as a strong support to the theory of special relativity. Now 
we have shown that this is once again a misinterpretation concerning about the classic physics. The 
transverse Doppler effect is definitely included in the exact solution (14) obtained using the principles of 
the classical physics. 


C. Field patterns 


The point source with harmonic charge is the simplest example for illustrating the radiation property. 
It can be treated as an imaginary source because it is not a constant and may violate the conventional 
charge conservation law. However, this is not fatal because we can combine two of them with opposite 
charges together to form a Hertzian dipole, the charge in which is conserved. Or we may simply consider 
it as a point source of scalar wave. Although the expressions for the wave velocity (13) and the frequency 
(14) are derived from a uniformly moving point source with harmonic time-varying charge, they can be 
used for predicting the propagation properties of the far fields of general sources in uniform motions. For 
general sources, the patterns of their fields may have much more complex directivities. However, the 
propagation velocity and the Doppler effect of their far fields are the same and can be respectively 
described by (13) and (14). 


As an example, we have derived the far fields of a moving Hertzian dipole. It is composed of two 
anti-phase harmonic charges with small distance / , as shown in Fig.3. The electric momentum is 
expressed by p,/cos(@,t,)p, where p is the polarization unit vector. The fields of the uniformly moving 


Hertzian dipole can be derived using exactly the same techniques in deriving the fields of the single 
harmonic charge. 


E(r,r) 


Po COS Ogh 


trajectory 
— pP, COS Of, 
Fig. 3. Trajectory of a moving Hertzian dipole. r is not depicted. 


We here consider the case that the dipole moves in a path perpendicular to its polarization direction, 
v-p = 0 , as shown in Fig.3. The scalar potential of the two charges is of the same form as (5), 


o.(r,t)=+ Pu 3 


dae Ril, py (eh) (19) 


where R,, n,, and ¢,, are respectively related to the charge +p, cos œt . The total scalar potential is 
then obtained by 


150 


olr) =n] AEDA (20) 


We let / —> 0 but keep p,/ constant. By applying the vector identities in deriving (7) [5], we get the 
final scalar potential for the moving Hertzian dipole, 


=> 


= : 7 COS —— = 
27e | R’(1--B) cyR(1-n-B) 
The vector potential for the Hertzian dipole is derived in a similar way, 


Pol JE EREDE osan) (â P)B+(1-û ‘B)p QD sin (@pt, | 


o(r,t) = pol (eee (a@t,)- us sausin(an)| (21) 


A(r,t) = 


47, 


c,R*(1—a-B) @R(1-a-By (al 


Note that the vector potential includes the contributions not only from the two moving charges but also 
from the moving dipole momentum. The fields can be derived from them rigorously with the same 


formulae as those in [4, 5]. Keeping the main terms that include the first order of (1/ R) , the far fields of 


the uniformly moving Hertzian dipole can be approximately expressed by 


nx(pxB)+nx(nxp 
E, (r,t) = Eo E- = 7 B) cos (œt ) (23) 
1. 
Bu, (r, t) = —nxE,, (r,t) (24) 
c 


0 


where E, = @ Pp! / (4z8,c> ) is a source-related constant. The propagation property is exactly the same 
as the uniformly moving single harmonic charge, except that there is a directivity factor related to the 
polarization of the dipole. For v =0, we can check that (23) turns into 


OE, 
By (yt) = iE 


(p—fF -p)cos(@,¢, ) 


which is the far electric field of a fixed Hertzian dipole. 


We will use three types of sources to illustrate the properties of the fields at 2 =0, 0.5, and 0.9. The 


electric fields for the time invariant charge, harmonic time varying charge are calculated with (7). The 
fields of the Hertzian dipole are calculated with (23). For motionless sources, the results are shown in 
Fig.4. The static charge simply generates a static field, as shown in Fig.4(a). For the sources with time 
varying charges, the wave velocity is c, in all directions, as shown in Fig.4(b) and (c). 


(a) (b) (c) 
Fig.4. The field patterns of the motionless sources. (a) Static charge. (b) Harmonic charge. (c) 
Hertzian dipole. 


The electric fields of the moving sources are shown in Fig.5 and Fig.6, in which red arrows indicate 
the moving directions. When the static charge moves, it generates a time-varying field in the space with 
a spectrum dependent on the moving velocity. Its field patterns contract along the moving direction. For 
moving sources with time varying charges, the anisotropic property of the wave velocity and the Doppler 
effect are clearly demonstrated in their fields. 


(a) (b) (c) 


Fig.5. The field patterns of the moving sources at 2 = 0.5. (a) Time invariant charge. (b) Harmonic 
time varying charge. (c) Hertzian dipole. 


(a) (b) (c) 


Fig.6. The field patterns of the moving sources at 2 =0.9 . (a) Time invariant charge. (b) 
Harmonic time varying charge. (c) Hertzian dipole. 

The radiation directivity of the Hertzian dipole is obviously different from that of the single moving 
harmonic charge. When the source velocity is close to c,, very strong radiations with very high 
frequencies can be observed in the moving direction of the Hertzian dipole, while radiations with lower 
frequencies can be observed in the opposite direction of the dipole. 

For £ = 0.9 , the wave velocity is shown in Fig.7(a). It is exactly c,,at 9=0,7, and is about 2.3c, 
at 0 =+/2. The normalized frequency is plotted in Fig. 7(b). In this case, 0, =132° and 228°. The 
red-shift area is much smaller than the blue-shift area. 
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Fig. 7. Moving Hertzian dipole for 6 =0.9 . (a) Wave velocity. (b) Normalized angular frequency. 


The normalized electric field is plotted in Fig.8(a). The radiation pattern is quite different from that 
of the motionless Hertzian dipole shown in Fig.8(b). Obviously, these results show that both the radiation 
pattern and the signal frequency of an antenna implemented on a fast-moving platform may suffer sever 
distortions and need to be carefully compensated. 
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(a) (b) 
Fig. 8. Radiation patterns of Hertzian dipole. (a) Moving at 2 =0.9 . (b) Motionless. 


Ill. Cherenkov Radiation and Electromagnetic Shock Waves 


As we can see, the fields we obtained are the exact solutions to the Maxwell’s equations. Since the 
Maxwell’s theory itself does not put any restrictions on the velocity of the electromagnetic waves, the 
formulae (7)-(14) are valid for all £, including £ >1. Therefore, we temporarily put aside the velocity 


limitation and analyze the properties of the causal solutions to the Maxwell’s equations when the sources 
move faster than light in the vacuum. However, we must observe the causality principle that the fields 
cannot be generated by the sources in the future. Explicitly, it is required that 4, should be real and ¢, <t. 


Without loss of generality, we choose the position of the source at time £ as the origin of the 
coordinate for the sake of simplicity, that is, we choose t=0 since x(¢) = vt. In order to guarantee that 
t, is real for t= 0, we derive from (10) that the range of the angle should satisfy 


B sin’ O<1 (25) 
which defines a biconical region. To meet the causality condition of t, < 0 , we check from (10) that 


cos0 <0 (26) 
Denote the critical angle as 
0, = z -sin™ (1/2) (27) 
It is obvious that the electromagnetic fields only exist in the conical region of 6, < 0 < z . In particular, 
at the edge of the conical zone, (1-û -B) — 0, the amplitude of the fields tends to become infinitely 
large for point sources. For general sources, the amplitude is finite but may be unusually large. The 
electromagnetic fields at the edge form a kind of electromagnetic barrier of very large amplitude. 


For a moving source with a constant charge, Cherenkov radiation occurs when £21 [10]. It is 


considered that Cherenkov radiation cannot occur in the vacuum since sources cannot move faster than 
the light in the vacuum. However, it can occur in a medium because sources can move faster than the 
light in that medium. Cherenkov radiation was discovered in 1934 in the experiment of bombarding the 
water with y -rays. Bluish light was clearly observed in the experiment. 


When the velocity of the sources with harmonic time varying charges increase to c,, then 2 =1, 


and 6, = 2/2, the wave fronts of the electromagnetic waves are in the same plane of the source. The 


source crosses the electromagnetic barrier [4, 10], similar to that of the sonic barrier that supersonic 
planes may encounter when they cross the sound speed. The patterns of the electric fields are shown in 
Fig.9. It is reasonable to predict that the electromagnetic shock waves may be observed in media. 


(a) (b) (c) 


Fig.9. Electromagnetic barriers generated by moving sources at 2 =1. (a) Time invariant 
charge. (b) Harmonic time varying charge. (c) Hertzian dipole. 


For superluminal sources, 8 >1, 0, > 2/2. The wave fronts lag behind the source. The fields are 


confined in the conical zone and cannot surpass the source. The solutions describe electromagnetic shock 
waves, just like the sonic shock waves. For point sources, at the edges of the shock waves, the fields tend 
to be infinitely large and form conical-shaped electromagnetic barriers. 


The electric field distributions for the uniformly moving point sources of the three types are shown 
in Fig.10. 


(a) (b) (©) 


Fig.10. Electromagnetic fields at 6 =1.5 . (a) Cherenkov radiation of a moving time- 


invariant charge. (b) Shock wave of a moving point source with harmonic charge. (c) Shock 
wave of a moving Hertzian dipole. 


At the electromagnetic barriers, not only the amplitudes of the fields are very large, the frequencies 
and the velocities of the fields also tend to become very large, as we can check that 


s(0.) B B cos 0, 


= o0 28 
-pP e J1- p° sin? 8, i ee 
c(0,)= -== > v (29) 


J1- p’ sin? 0 
The superluminal sources radiate high energy rays with velocity much larger than c, at the 
electromagnetic barriers. The observer will experience an electrmagnetic boom when he is swept by the 
electromagnetic barriers of the superluminal sources. However, for an observer outside the shock wave 


zone of a superluminal source, the superluminal source is invisible to the observer. We may consider the 
superluminal source as an electromagnetic blackhole to the observer. 


The wave velocity, the normalized frequency, and the amplitude of the fields in the shock wave 
region in the far region for 2 =0.9 are plotted in Fig.11. As can be seen, the fields in the middle area of 


the shock wave conical region have the largest red-shift and propagate with velocities slightly larger than 
Gye 


8 [°] 
Fig. 11. The velocity, the normalized frequency, and the normalized field amplitude in the shock 
wave region for B=0.9. 


IV. Conclusions 


It is often stated in textbooks that, according to the rules of the classic physics, the wave velocity 
and the velocity of the sources are added like vectors; the light velocity in the vacuum is constant in all 
directions and is independent upon the light sources; the transverse Doppler shift is a pure relativity effect. 
These are three longstanding misinterpretations concerning about the classic physics. The exact solutions 
of the electromagnetic fields of uniformly moving point sources in the vacuum reveal that the far fields 
are spherical waves with velocity depending on the propagation direction. The wave velocity and the 


moving velocity of the source satisfies c = c, i V1-{’ sin’ 0 . This relationship clearly demonstrates that 


the two velocities cannot be added like vectors under the principles of the classic physics. Moreover, 
transverse Doppler effect is intrinsically included in the solutions. 


Although superluminal sources have not been officially confirmed [11, 12, 13], it does not discourage 
us from performing mathematical analysis on the causal solutions of the fields of sources moving faster 
than the light in the vacuum. The electromagnetic shock waves are naturally introduced according to the 
solutions, together with a reasonable hypothesis that the superluminal sources can be considered as 
electromagnetic blackholes to observers staying in regions outside the shock wave zone. 
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